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We investigate the thermodynamic properties of the frustrated bilayer quantum Heisenberg an- 
tiferromagnet at low temperatures in the vicinity of the saturation magnetic field. The low-energy 
degrees of freedom of the spin model are mapped onto a hard-square gas on a square lattice. We 
use exact diagonalization data for finite spin systems to check the validity of such a description. 
Using a classical Monte Carlo method we give a quantitative description of the thermodynamics of 
the spin model at low temperatures around the saturation field. The main peculiarity of the consid- 
ered two-dimensional Heisenberg antiferromagnet is related to a phase transition of the hard-square 
model on the square lattice, which belongs to the two-dimensional Ising model universality class. It 
manifests itself in a logarithmic (low-)temperature singularity of the specific heat of the spin system 
observed for magnetic fields just below the saturation field. 
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I. INTRODUCTION 



The quantum Heisenberg antiferromagnet (HAFM) on 
geometrically frustrated lattices has attracted much at- 
tention during last yearsJ^ Besides intriguing quan- 
tum ground-state phases at zero magnetic field those 
systems often show unconventional properties in finite 
magnetic fields like plateaus and jumps in the magne- 
tization curve, see e.g. Ref. y. The recent finding 
that a wide class of geometrically frustrated quantum 
spin antiferromagnets (including the kagome, checker- 
board and pyrochlore lattices) has quite simple ground 
states in the vicinity of the saturation field)^ namely in- 
dependent localized-magnon states, has further stimu- 
lated studies of the corresponding frustrated quantum 
antiferromagnets at high magnetic fields iSiSiL2iS In par- 
ticular, the low-temperature high-field thermodynamics 
of various one- and two-dimensional frustrated quan- 
tum antiferromagnets which support localized-magnon 
states, can be discussed from a quite universal point 
of view by mapping the low-energy degrees of freedom 
of the quantum HAFM onto lattice gases of hard-core 
objects iSiSiSiiS For instance, the kagome (checkerboard) 
HAFM in the vicinity of the saturation field can be 
mapped onto a gas of hard hexagons (squares) on a tri- 
angular (square) lattice The exactly soluble hard- 
hexagon model exhibits an order-disorder second-order 
phase transitioniii The hard-core lattice-gas model cor- 
responding to the checkerboard HAFM consists of large 
hard squares on the square lattice with edge vectors 
fli = (2, 0) and a2 = (0, 2) (i.e. there is a nearest- 
neighbor and next-nearest-neighbor exclusion). For the 
latter model no exact solution is available, but most likely 
there is also an order-disorder phase transitioufi^ The ex- 
istence of a phase transition in the hard-hexagon (large- 



hard-square) model would imply a corresponding finite- 
temperature transition of the corresponding spin model 
near saturation provided the low-temperature physics is 
correctly described by the hard-core lattice-gas model. 
However, at the present state of the investigations no 
conclusive statements for the kagome and checkerboard 
antiferromagnets are available, since both models ad- 
mit additional degenerate eigenstates not described by 
the hard- hexagon/large-hard-square model A?-^- Further- 
more for these spin models precise statements on the gap 
between the localized-magnon ground states and the ex- 
citations are not available. Therefore, the effect of addi- 
tional ground states and the excited states on the low- 
temperature thermodynamics remains unclear. 

The motivation for the present paper is to find 
and discuss another two-dimensional frustrated quantum 
HAFM, for which a hard-core lattice gas completely cov- 
ers all low-energy states of the spin model in the vicinity 
of the saturation field and where all excitations are sep- 
arated by a finite energy gap. For such a spin model 
one can expect that an order-disorder phase transition 
inherent in the hard-core lattice-gas model can be ob- 
served as a finite-temperature phase transition in the spin 
model. It might be worth to note that such a phase tran- 
sition of course does not contradict the Mermin- Wagner 
theoremi^ that forbids magnetic long-range order (break- 
ing the rotational symmetry) for the two-dimensional 
Heisenberg model at any non-zero temperature and at 
zero field. 

A spin model which satisfies these requirements is a 
frustrated bilayer quantum HAFM. The investigation of 
the bilayer quantum HAFM was initially motivated by bi- 
layer high- Tc superconductorsi^ and has been continued 
till present time, see e.g. Ref. ^3 and references therein. 
Below we will illustrate that the corresponding hard-core 
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spins of length s reads 



FIG. 1: (Color online) The frustrated bilayer antiferromagnet 
(black lines). The vertical bonds have the strength J2 whereas 
all other bonds have the strength Ji . The trapping cells (ver- 
tical bonds) occupied by localized magnons are shown by fat 
lines. The auxiliary square lattice is a simple square lattice 
filled by hard squares (indicated as red squares) which corre- 
spond to localized magnons. 



lattice-gas model is a model of hard squares on a square 
lattice, however, in difference to the checkerboard lattice 
with smaller hard squares with edge vectors ai — (1,1) 
and 02 — (—1,1) (i.e. there is a nearest-neighbor exclu- 
sion, only, cf. Fig. This hard-square model exhibits 
an order-disorder phase transitioniiiiiSiiS In the context 
of different universality classes discussed in Ref. [ol the 
frustrated bilayer quantum HAFM is the first example of 
a spin system which belongs to the universality class of 
(small) hard squares. 

The paper is organized as follows. First, we specify the 
frustrated bilayer model and illustrated the correspond- 
ing localized- magnon states (Secnj. Then in Sec. lIIII we 
calculate the contribution of the independent localized- 
magnon states to the thermodynamic quantities using the 
hard-square model. We compare our results with exact 
diagonalization data for finite spin-1/2 Heisenberg sys- 
tems of up to TV = 32 sites. Finally, in Sec. II VI we report 
the low-temperature high-field thermodynamic quanti- 
ties obtained on the basis of Monte Carlo simulations for 
hard squares focusing on the (low-)tempcrature depen- 
dence of the specific heat in the vicinity of the saturation 
field. 



II. THE FRUSTRATED BILAYER 
ANTIFERROMAGNET AND INDEPENDENT 
LOCALIZED MAGNONS 



To be specific, we consider the nearest-neighbor 
Heisenberg antiferromagnet in an external magnetic field 
on the lattice shown in Fig.^ This lattice may be viewed 
as a two-dimensional version of the frustrated two-leg lad- 
der considered in Refs. 20 21 22 (for some similar models 
see Ref. l23j) . The Heisenberg Hamiltonian of N quantum 
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Here the sum runs over the bonds (edges) which connect 
the neighboring sites (vertices) on the spin lattice shown 
in Fig. ^ Jnm > are the antiferromagnetic exchange 
constants between the sites n and m which take two val- 
ues, namely, J2 for the vertical bonds and Ji for all other 
bonds. A > is the exchange interaction anisotropy pa- 
rameter, h is the external magnetic field, and — 
is the 2-component of the total spin. In our exact diago- 
nalization studies reported below we will focus on s = 1/2 
and A = 1. 

We note that the Hamiltonian commutes with the 
operator and hence we may consider the subspaces 
of its eigenstates with different values of separately. 
Evidently, the fully polarized state |0) = |s, . . . , s) is the 
eigenstate of the Hamiltonian with = Ns and can 
be considered as the vacuum state with respect to the 
number of excited magnons. This state is the ground 
state for high magnetic fields. 

Consider next the one-magnon subspace with ~ 
Ns ~ 1. The one-particle energy is given by 
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(1) 



-s(J2 + A(8Ji-f J2)) -I- /i and A 



(2) 



s (4 Ji (cos + cos ky) + J2 

ka — 2Trna/Ma, Ua = 1,2, 



A(8Ji + J2)) + h (here 
.,M„, x,y, M^My = 

N/2). Obviously, the excitation branch A^^-* is disper- 
sionless and it becomes the lower one when J2 > 4Ji. 
Throughout this paper we assume J2 > 4Ji. Then the 
saturation field is given hy hi = s ( J2 + A (8 Ji + J2)). 

The A'^/2 dispersionless one-magnon excitations can 
be written as localized excitations on the N/2 vertical 
bonds, i.e. |1) = |lm)„|s, . . . , s)e is an eigenstate of ^ in 
the subspace with S"^ = Ns — 1 with the zero-field eigen- 
value EpM - ei, where Epu iNAs^Ji + NAs'^J2/2 
and ei = s ( J2 + A (8J1 -I- J2)). In |1) the first part is 
the localized one-magnon excitation on the vertical bond 
number v, i.e. |lm)^ — 2^^!"^ (|s, s — 1) — |s — 1, s))^ and 
the second part |s, . . . , s)e is the fully polarized environ- 
ment. 

We pass to the subspaces with — Ns — 2, Ns — 
3,...,Ns — n„iax, where nmax = A^/4. We can easily 
construct many-particle states in these subspaces using 
the localized- magnon states. Explicitly the wave function 
of n independent localized magnons has the form 



\n) = |lm)„i . . . |lm),„^ . . . \lm)y^ \i 



(2) 



It is important to note that any two vertical bonds Vi 
and Vj in Eq. |(2Jl where localized magnons live are not 
allowed to be direct neighbors. The energy of the n inde- 
pendent localized-magnon state Q in zero field ft, = is 
En = EpM — nei. Since n independent localized magnons 
can be put on the bilayer in many ways the eigenstates 
^ are highly degenerate. We denote this degeneracy by 
gN^n), that is the number of ways to put n hard squares 
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on a lattice of = N/2 sites (see Fig.^. According to 
Refs. the independent localized-magnon states Q 
are the states with the lowest energy in the correspond- 
ing subspaces. Moreover, they are linearly independent 
(orthogonal type in the nomenclature of Ref. [ij) and 
form an orthogonal basis in each subspace^ Due to their 
linear independence they all contribute to the partition 
function of the spin system. 

In the presence of an external field the eigenstates (0) 
have the energy En{h) = E'fm — hsN — n{ei — h). At 
the saturation field, h = hi = ti, all they are ground 
states and the ground-state energy £'„(/ii) does not de- 
pend on n. As a result the ground-state magnetization 
curve exhibits a jump at the saturation field. This jump 
is accompanied by a preceding wide plateau, where the 
width of this plateau can be obtained following the ar- 
guments given in Ref. 22 and from finite-size data. We 
find for s = 1/2, A = 1 a plateau width of /ii — /12 = 4Ji. 
This plateau belongs to the two-fold degenerate ground 
state with maximum density rimax = N /A of localized 
magnons, the so-called magnon crystal,** where all local- 
ized magnons occupy only one of the two sublattices of 
the underlying square lattice. 

Furthermore, the degeneracy of independent localized- 
magnon states at the saturation field W = X^n^cT 9N{n) 
grows exponentially with the system size N that im- 
plies a nonzero ground-state residual entropy S = 
fclimAr^oo(lnW/iV) = 0.2037. . .fc (see Ref. Eland also 
below). Due to their high degeneracy the independent 
localized-magnon states are also dominating the ther- 
modynamic properties at low temperatures for magnetic 
fields around the saturation field as we will dicuss in de- 
tail below. 



III. HARD-SQUARE MODEL 

We want to calculate the contribution of the indepen- 
dent localized magnons to the canonical partition func- 
tion of the spin system. 



Zira{T,h,N) = ^ 5Ar(n)exp 



Tl = 



kT 



exp 



Eym — hsN 
kf 



^ gAr(n)exp (^^n] , (3) 



where fi = ei — h = hi — h. It is apparent that 
(7Ar(n) is the canonical partition function Z{n,J\f) of n 
hard squares on a square lattice of = A^/2 sites, 
whereas E{T,fi,N') = Yl^^a 9N{n)ex.Y){^n/kT) is the 
grand canonical partition function of hard squares on 
a square lattice oi N = N/2 sites and fi is the chemi- 
cal potential of the hard squares. As a result we arrive 
at the basic relation between the localized-magnon con- 
tribution to the canonical partition function of the spin 
model and the grand canonical partition function of hard- 



square model, 

ZUT,h,N) = exp (^-:^^M_^) S(T,M,A/-). (4) 

Eq. Q yields the Helmholtz free energy of 
the spin system Fu„(T,/i,7V) = -kT\wZx^{T,h,N), 
whereas the entropy S*, the specific heat C and 
the magnetization M ~ (5^) are given by the 
usual formulas, Sin,{T,h,N) = -dFi„,{T,h,N)/dT, 
Cin.iT,h,N) = TdSUT,h,N)/dT, MUT,h,N) = 
sN - kTd In s(r, ^l, TV) /^^l. 

We use exact diagonalization data to check this picture 
for s = 1/2 isotropic Heisenberg system ^ oi N = 1& 
and = 20 sites (full diagonalization) and A^ = 32 (only 
in the subspaces with — 16, . . . , 11) imposing periodic 
boundary conditions. We fix the energy scale by putting 
Ji — 1. For the vertical exchange bonds we consider 
J2 = 4, 5, 10. 

First we compare the degeneracies gN{n) of the local- 
ized n-magnon states calculated for spin systems of sizes 
A^ = 16,20,32 with the corresponding values Z{n,N/2) 
of the hard-square model. As expected we find com- 
plete agreement between both models for J2 > 4Ji. 
As an example we give here the numbers for A^ = 20: 
g2o(n) = 1, 10, 25, 20, 10, 2 for n = 0, 1, 2, 3, 4, 5. 
For J2 = 4Ji the spin system has one extra state for 
n = 1, i.e. in the one-magnon sector, since the dispersive 
mode K'^^ and the dispersionless mode A^^-* are degener- 
ate at k = (tt, tt). 

To estimate the relevance of excited states, not de- 
scribed by the localized-magnon scenario, we have de- 
termined the thermodynamically relevant energy sepa- 
ration Ados between the localized-magnon states and 
the other states of the spin system by calculating the 
integrated low-energy density of states at saturation 
field. We define Ados as that energy value above the 
localized-magnon ground-state energy, where the contri- 
bution of the higher-energy states to the integrated den- 
sity of states becomes as large as the contribution of the 
localized-magnon states. For both values J2 = 5 and 
J2 — 10 we find Ados ~ 1 independent of the size of the 
system A^. We can expect that the contribution of the 
localized-magnon states to the partition function is dom- 
inating for temperatures kT significantly smaller than 
Ados- 

In Fig. 121 we present our results for the magnetization, 
the specific heat and the entropy (panels from top to 
bottom) for the spin system of size A^ = 20 (triangles 
and diamonds). We compare those data with the corre- 
sponding data taking into account only the contribution 
of independent localized-magnon states, described by the 
finite hard-square model of M = 10 sites (thin dotted 
lines). Up to kT w 0.1 both data sets coincide, demon- 
strating that the hard-square description perfectly works 
at low temperatures. But we observe good agreement 
also for higher temperatures up to kT k, 0.3, indicating 
that the independent localized magnons still dominate 
the thermodynamic quantities. Further increasing kT, 



4 



the high-energy states more and more contribute to the 
partition function and the hard-square description loses 
its vahdity. Note that an identical statement can be made 
for N = 16 {J\f ^ 8). 




-1 1 
(h - hi)/(kT) 



FIG. 2: The magnetization {sN~M{T, h, N))/N (a), the spe- 
cific heat C{T,h,N)/kN (b) and the entropy S{T,h,N)/kN 
(c) as functions of {h — hi)/kT. Symbols correspond to the 
exact diagonalization data (A'' = 20, J2 — 5, kT — 0.2 (tri- 
angles), kT = 0.3 (diamonds)). Thin lines correspond to the 
results for a small hard-square system A/" = 10 (thin dotted 
lines) and TV = 64 (thin solid lines). The exact diagonaliza- 
tion data for the temperatures up to kT = 0.1 coincide with 
the corresponding data for the TV = 10 hard-square system. 
We also show by thick lines the Monte Carlo simulation results 
for the magnetization and the specific heat for a hard-square 
system of sizes up to 800 x 800 sites. 



IV. LOW-TEMPERATURE 
THERMODYNAMICS AROUND THE 
SATURATION FIELD 

Now we discuss the low-temperature high-field thermo- 
dynamics of the frustrated bilayer quantum HAFM using 
the results for the hard-square model. From Fig. it 
is obvious, that fixing the temperature to a sufficiently 
low value the spin system can be driven through a phase 
transition by increasing the magnetic field h towards the 
saturation field hi. On the other hand, we can fix the 
magnetic field slightly below the saturation field and vary 
the temperature. Then the phase transition is driven by 
the temperature and the specific heat exhibits a singu- 
larity at a critical temperature Tc{h), see Fig. |21 where 
we show the results of a Monte Carlo simulation for large 
hard-square systems (periodic cells with up to 800 x 800 
sites and 3-10^ steps) for the specific heat for two val- 
ues of the magnetic field. The data clearly indicate a 
phase transition which occurs in the hard-square model 
at Zc = 3.7962. .., i.e. at {{hi-h)/kT)c = Inz^ « 1.3340 
which yields kTc{h) f» {hi — ft,)/1.3340. The correspond- 
ing order parameter is the difference of the density (of 
hard squares or localized magnons in hard-square or spin 
language, respectively) on the A- and B-sublattices of the 
underlying square latticeiii For {{hi — h)/kT)c < Inzc 
(i.e. for T > Tc{h), h < hi) both sublattices are 
equally occupied, but for {{hi — h)/kT)c > Inzc (i.e. for 
T < Tc{h), h < hi) one of two sublattices is more occu- 
pied than the other. Therefore in the ordered phase the 
translational symmetry of the spin (hard-square) system 
is broken. Finally, at T = only one sublattice is occu- 
pied and the other is empty and the ground state of the 
spin system is a magnon-crystal state, see Sec. 

We can estimate the critical temperature for a fixed 
deviation of the field from the saturation value using the 
above given expression for kTc{h). For \ — h/hi = 0.02 
we find kTc/hi « 0.0150, whereas for 1 - h/hi = 0.01 
we have kTc/hi w 0.0075 and for the spin system 
with J2 = 5 we find for h = 8.91Ji {h = 8.82Ji) 
kTc « 0.0675Ji {kTc « 0.1349Ji). Such temperatures 
are within a temperature range where the hard-square 
description for finite systems works perfectly well (see 
Figs. 13 EJ- However, one may expect that the scenario 
of the phase transition may hold also at temperatures, 
for which the localized-magnon states are still dominant 
but also higher-energy states of the spin system not de- 
scribed by the hard-square model contribute to the par- 
tition function. 

We mention that the hard-square model belongs to 
the two-dimensional Ising model universality clas3ii*2S*S& 
with the critical exponents /3 = 1/8 for the order param- 
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J2=5,N=20, 1.02hi 
J2=5,N=20, l.Olhj 
HS, 10 sites, 1.02hi 
HS, 10 sites, l.Olhj 
HS MC, 1.02hi 
HS MC, l.Olhi 
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eter and a = for the specific iieat, i.e. tiie specific tieat 
shows a logarithmic singularity at the critical point. Note 
that this universality class is different from the one of the 
hard-hexagon modelpii Thus, the low-temperature peak 
(singularity) in the temperature dependence of the spe- 
cific heat in the vicinity of the saturation field is a spec- 
tacular sign of highly degenerate independent localized- 
magnon states of the frustrated bilayer quantum HAFM. 
Their ordering leads to the hard-square type peculiarity 
just below the saturation field. 

Another thermodynamic quantity of interest is the 
curve of constant entropy as a function of magnetic field 
and temperature. Since hard-square description implies 
the dependence of the entropy only on (hi — h)/kT this 
curve is similar to an ideal paramagnet. Therefore, the 
considered spin system is expected to exhibit a large 
magnetocaloric effect in the vicinity of the saturation 
field&La 

Finally we note, that the effect of the localized 
magnons is a pure quantum effect which disappears as 
s increasing approaches the classical limit s ^ oo^ 



FIG. 3: The temperature dependence of the specific heat 
C{T, h, N) /kN for h = l.Ol/ii (broken lines, diamonds) and 
h = 1.02^1 (solid lines, triangles) (upper panel) and for 
h = O.99/11 (broken lines, diamonds) and h = 0.98/ii (solid 
lines, triangles) (lower panel). The results of Monte Carlo 
simulations for a hard-square system of sizes up to 800 x 800 
sites are shown by lines. We also report the corresponding re- 
sults for the hard-square system with — 10 (open symbols) 
and the exact diagonalization data for the finite spin system 
with AT = 20, J2 = 5 (filled symbols). 
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